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ABSTRACT 


Intensity scintillation has been studied for a He-Ne 
(6328A) laser beam propagating over a 4.05 kilometer hori- 
Zoncal path tm the marane boundary layer. The results 
fall in three areas. (1) The covariance of the logarithm 
of the laser beam amplitude increases significantly with 
@eereasing detector aperture diameter for apertures less 
tie hemita@menesielma zone Size. this effect appears to 
be dependent on the magnitude of the refractive index 
Senucture constant, Cy (2) The dependence of Cy on 
meine) FOr suUnseab Lemmecreerological eenditions was 
verified using the laser beam in determining Cy. (3) Fre- 
G@uency analysis performed on scintillation data indicates 
Mmiapethe increase of power in higher frequencies is more 
dependent on the wind velocity perpendicular to the propaga- 
Miomepath than on detector aperture diameter for apertures 


fier chan the £irst Fresnel zone size. 
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I. INTRODUCTION 


The structure of the atmosphere is known to modify a 
feacer Deam propaeatineg through if. The light may be 
Posorped dic to atiospmerie gases and pollutants, scattered 
by molecules, electrons or particles, or undergo refraction 
due to turbulence caused by wind and temperature and pres- 
SUncme raGdmemicwmem erraction €fLiects can be further divided 
into beam wander, beam spread and beam scintillation. 
Oyerall, these effects comprise a significant source of 
errors and distortions in laser communications, radars and 


Comtrol systems. 


wee ROJERCIMAIS TORY 

Pe loviman Titereisceiplimary ¢roup of professors at the 
fewal Postgraduate School formed an Electro-Optics group for 
mie purpose Of pursuing advanced research in a wide variety 
eee clectro-optic areas, ~ihis grotsp has members from Physics, 
Meteorology, Oceanography, Mechanical Engineering, Electrical 
Engineering, Ger oliteterleanoimecring, A team has been 
investigating aspects of laser beam propagation over the 
Seed surtace under a Navy research contract. Previous 
theses by physics students have dealt with investigation of 
meme Modulation transfer function [Ref. 8] and with verifica- 
tion of the Kolmogoroff turbulence model and comparison of 
Cy, as determined by optical methods to that determined by 


micro-meteorological methods [Ref. 6]. This thesis 






undertakes the investigation of laser beam scintillation in 
tHemeCentext of averture-averaging effects, height dependence 
and frequency spectra. Other work being done simultaneously 
involves the determination of laser beam wander and beam 
spread [Ref. 9]. The work thus far has been with 6328 


Angstrom He-Ne lasers. 


B. SCINTILLATION 
ieee ese: ptron 

Saunt) | lationesmecancmetwinkling'' or, sparkling" and 
momrcadily observable if one looks at distant astronomical 
pocdtes., ithe twinkling 1s due to changes in the refractive 
midex of the atmosphere. 

A simple model of laser beam scintillation can be 
described by considering the propagation path to be made 
Meeot small parcels of unauform density (and, therefore, 
Selscaie rertactime  lndex) see lent ray is bent slightly 
Meaeeach parcel» dependane on the value of the parcel's 
refractive index, n. Meo varlation in path between 
adjacent wave fronts causes interference patterns to be 
formed. The effect at a point receiver is a time varying 
mmECNSIty signal. 

fe LieOrTy 

tie winictenteranaonness of turbulence leads to a 
aeoctotiealemodeling Gt scintillation. The refractive index 
variations along the propagation path modulate the intensity 


of a beam in a multiplicative manner; thus they modulate the 






logarithm of the intensity (and amplitude) in an additive 
Manner. the variation of the logarithm of the amplitude is, 
therefore, the sum of the random perturbations induced along 
Pipa essa Consequence Of the central-limit theoren, 
the variations of the log-amplitude should follow a normal 
frstribution. 

ics io muiesmOnmm tne stUrbDuhent mixing of the 
atmosphere are provided by Kolmogoroff's analysis [Ref. 1]. 
ao ieon sOretinnuUlemec qismaws Che mean sqWare difference 
of the refractive index, measured at two points separated 


Be a2 distance d, to be written as 


i 


(1) p,(4) = <{n(r) - n(r+a)]7> 


i 
Qu 


D, (4) 


D Cd) Ewe welCdeEne NetmIcLive wnaex structure function. 
E>. : , 

Cy Poetics cihdechivec Inge>wstnhUeture Constant and param- 

feerizes the total amount of energy in the turbulence. 


Tatarski [{Ref. 2] has formulated the three-dimensional 


eeectral density of the turbulence as 


(2) $6 (K) = 0.033 ¢ ay-1i/s 


ee = 2 
N exp[-K /K_ ] 


Here K = 2m/& and is the spatial wave number corresponding 
meee, the coherence length of a uniform turbulent parcel of 
aay . Kn is the spatial wave number corresponding to the 


inner scale of turbulence, Le. This length is the lower 





MikmtmOn tice inertial sub-range of turbulence. For values 
of & less than Ros Viscous dissipation of energy begins to 
aoniiate., [he outer scale of the inertial sub-range is L 
miascorresponds to the range in which energy is added to 
ere turoulence Spectrum, The Kolmogoroff theory applies 
Smiy to the region &o<k<hi. iwervaiue of de has been 
weimaced tO be om tire order of a few millimeters while Oe 
mmole Order Of Meters. 


thewsedlamwewave equation 1S given by 


(3) pon @ once St 


assuming that any polarization fluctuations can be neglected, 
and that \<<h A is the wavelength of the propagating beam 
mack tTlow 1S Cqual to 2Z7/A. Using Rytov's method [Ref. 3} 


pod letting E = exp[W], (3) becomes 


(4) em @ewet nica’ = 0 


Meplying a first order perturbation of 


<a 
HI 
<— 
=) 
a 
— 
— 


and 


to equation (4) leads to 


2 2 _ 
eS) Vv, + 2VW_*V, + 2k iy = tt 





Mire Solution 15 
(6) b. (2) = 1/E,(r) f 2k’n, (%,)E,(£')G(r-4') dV" 
1 Dee. Tor 6 


loco aecevitimecntcrmone cy ochmeltzer [Ref. 4] 
obtained the statistical expression for the log-amplitude 
Sovariance, cf (a), eter rt edmnet., o| applies to a 


spherical wave propagating in a random medium with the result 


WE co 
Mac” (d) = k-/en Rei f ds fF 6 (K*’*) J. (K*/2as/z) 
aN ' 0 n 0 
he ieeeercp (ks(z-s)/iKs]dk 

and 

2 7 a ae mor ll) 6. 2 
es ) Co, (d=0) = COA = 0.124k Zz Cy 
where 

e api (Aye in <A> Bs 
QA n n ] 


Pim JO Sit TOmm meon( Sy) etne assumption Nees 
is made allowing the binomial expansion of the exponential 
and the dropping of higher order terms. This is a weak 
restriction at optical wavelengths. 


The Log-intensity covariance is then 


(9) Oo Ene ee 





waren, leads to 


7/12,-11/12 C 


(10) G. Si, ae 1 sf 


N 


The following is a summary of the assumptions and 
Meetrictions that apply toe this result: 
fe ihe turbulence energy spectrum can be described 
by the empirical Kolmogoroff model (equation 1). 
Meal Oolirizati One mine tlatrons may be neglected. 
Meomly the turbulent eddies in the inertial sub- 
MaOMaomeOMtriIbueeesienii1cantly to scintillation. 
Tole Wavctenathmemmonagarecdgisemuch less than the 
Tone Gb Sed heOROL eto lence, hoe 
TiC mEutCtUdatiomemim che refractive index are 
much smaller than the mean value. 
Ox n?z/a *<<i 
ee lurpudence 1S 9rsotrepic and homogeneous. 
ee line wave propagated 1s spherical. 
pp (Ne wcletector ob fhe imGensity changes is a point 
_ detector located on the beam axis. 
mi ihe first order perturbation technique used in 
the Rytov method adequately models the wave 
propagation. 


mee ine beam is propagated over a horizontal path. 


10 





it. PROBLEM 


fe APERTURE AVERAGING 
Micepbeste .olmegorerr model depends on the experimental 
aeeeetor beime a pOant reeeirver. The general effect of having 
MeetireTGCtECCEPON SUGLaee 1S to reduce the degree of inten- 
moeriWettattons sensed, to reduce the frequency of observed 
changes in intensity and to increase the average signal 
Sexength. 
ie meOory 
For atmospheric laser propagation, the amplitude 
meirelation distance 1S of the order of Cua tne ia is & 
Fresnel zone size [Ref. 2]. 
tenet teornote 7 Oncmsmacimmed on the collection 
feemuulre as Chateradius R at wWiwen the geometric path from 
men sOurce to R and from R®to the detector 1s one-half a 
Peer lonerh tonger tham the axzal path. if the diameter of 


Ve 


fame recelving aperture is greater than Oa tee Wan 
Sentcain wave front sections with fluctuations of opposite 
Sign so that. the overall light flux will fluctuate weakly 
Ssempared to that seen Pyeetreaperture with a diameter iess 
oman the first Fresnel zone size. 

a. Wave-Optics Approach 


Pome sem incelgeht intensity, The total 


imeeict flux through the aperture is 


Vide 





(11) P = sf I(x,y) dxdy . 
A 


fem itictuetrens im P are defined as P' = P - <P> , and 
(el 2) Pee iel oe) axdy 

A 
where I' = [ - <{> 


The mean square fluctuations are 


Z 


Mey <P' "> = Ssssf <I'(x,,y,)1' @,,y,) 7 dx, dx,dy, dy, 
ALA 
eZ 
Let Wy = Cie aimue22" 2) = W, (rv) 


where r is the distance between the two points. 
Momence (xX ,y) SUCH tENatwite's zero OUtSide the aperture and 


irom tts surface. Also define 


feawation (13) may then be written as 


Cc 


7 oo ; 
meq) < PI > — if Wy Cao aye PS a INES ig BSG 


Now let 


Oo 


(15) H(x,y) eet 8 Y) x, ty, 


dRZ 





Peeiating Hofor a Circular aperture of radius R 


: Zz 2 
F(x) 54) = 1 for ie? ee 


: = : 2 - 2 Zi 
F(x-x,,y-y4) splits: ong (x x4) * Cr Yy? Ss R 
at TC OOnmise a cluche Of TFadius R centered at the 
Panel ine second ts ascirecle of radius R centered at 
(x45) Ott om on clcmaiLcanwonunica by the intersection 


Memthe two circles. 


mio) H(x,y) = H(r) IR? farccos (1/2R)-1/2R(1-r2/4R-) 1/2} <2R 


=) ee ea 


Substituting equation (16) into (14), 


2R 
(17) <P'*> = 20 f W,(r) H(r) r ar 
0 


H(r) equals 1/2 i Oraaa VMeCpOlmemapenrture), and goes to 


mewomtor © = 2R. tJatarski [Ref. 2] states that for 


el 2 


R<< (AZ) Rite Shine eakom W > @) Pmaopmoximately equal to 


W, (0) ewok eerOretiis case 


(18) <Pi4> = (aR2) %<1-i>? 


Wetarsion a2 ls omisesce thee rat io 


13 





2 





<P17> 
where Oh) = 7 
<p> 
. 
and mes ae ee 
OG as Presae) . 
alee 


iivesratio G(R) shows by what factor the fluctuations of the 
eon flux through an aperture of radius R are less than 


mmose through a vanishingly small aperture. 


2 


ZR 
(19) G(R) = 4/7R” Sf W, (1) /W; (0) [arccos (r/2R) 
0 


52 


EPPyOR lor 74n-)-’ “trdr 


G(R) equals unity for R equal to zero and monotonically 
decreases as R increases. 

omeeciitarciomevammate (17), the functions W,(@) 
and W, (0) must be known. Fried [Ref. 10] has written G(R) 


fons 6) in terms of aes picmeoveanuance Of the log-amplitude. 
2 . 
Cae seb aey  o 2-1 (*(X,,y,) - <2 J> 


Where £(x,y)' is the log-amplitude and is related to l1(x,y), 


Mme intensity, by 


I(x,y) = 1, exp [(22(x,y)] 


As noted above, 2&(x,y) is Gaussian random variable. Using 
his data to get values for CrINGE) amarassuming diitferent 
values for ci. (0), Porecdmerlemlaccd G(k) for normalized 


aperture diameters D/(42/k) 1/2. 


14 





mot wenn Oor tte umermalized aperture diameter, 
Mespnedicts G@y showlld vary as the anverse square of the 


aperture diameter. 


a2 ) G(R) = p2 Mecoustant. 


Miers ki | Rete) Zi), On the other hand, predicts 


(21) CaeenDy 2az)y. 7) (/° 


mot LD >> 2 On and small values of Coy (see equations (8) 
oad (9)). 


iMicmiimem~cerrTCeNON Ol, £n1lsS theory 1S that it 


mepliecs tO intinite plane waves for horizontal propagation 


paths. 
bemoconerrteal Optics Approach 
Young, an astronomer, approached the problem 
Heep geometric Optics [Ref . 11]. For an aperture much 


merci tian che inner Scale ot turbulence or the diffraction 
Beale, Oe. feolelemermsomercs as valid. Assuming the 
meeature is small compared to the outer scale of turbulence, 
mie hKOlmoporoff isotropic turbulence spectrum produces a 
shadow pattern with modulation power density proportional 

to eK, where K is the spatial wavenumber. The fractional 
modulation power of K remaining after integration over a 


Bebreuler aperture of radius R is 


£(RK) = [2J, (RK)/RK]? 


15 





Mhere £(RK) 1S the aperture filter function. 


ihemtetal modulation power then is 


ie?) p . sr 22K /3e(RK) K aK dé 
- anaep 2 ¢ K72/3 J, ° (RK) dk 
0 
Letting U = RK gives 
3) _ 25R/3 y 2/5 J, cu)au = 66 von es 


Inspection of the integrand of equation (22) shows that it 
meels Olt tTapidly for large K; most of the contribution comes 


mom a Narrow range of kK ~ 1/R. Observations of stellar 


femnitcillation contirm P ~ Ps iiwete. largemapentine 
range. 

oupeeoiivereativemor this result is for 
R= (az)i/? 


Poe EWG gee me 11/6 


the result for Hits OoiewaMmiced turbulence (i1.e., very 
Small apertures). This is unexpected for a theory fornu- 
meoed On geometric optics. 
Pie ocotmilecrtionvon this theory is that R must 
be >>(az)l/2, 
Pew Miteintea Experimental Data 
Pit odweMevers amawheilster made laser beam 


scintillation measurements with 21 collection apertures 


iG 





Peo nomenon iimerer to One meter in diameter [Ref. 7]. 
These measurements were made with a 6328A He-Ne laser over a 
mloping 8 Km patm. The terrain was rolling so that the beam 
jedmeam average ne1peht above the surface of greater than 15 
Meters. Data accumulation periods were one minute long with 
a time separation between data periods of four minutes. 
Peeure I is a reconstruction of their data. oe decreased 
with aperture diameter from 0.438 at 1 mm to about 0.378 at 
hoomm aperture, or about 86 percent of peg (1 mm). From 


fem to 100 mm the drop was from 0.382 to 0.220 such that 


Z A 
LA LA 


mooted’ Ss theory as presented in Ref. 10 and Section A.l.a. 


C CLeCeim)) = W557 ¢ (10 mm). The data did not dispute 
above. His data should be evaluated in terms of possible 
sources of error in Cre measurement due to non-simultaneous 
receiver background compensation, and the assumption that 
cr was the same over the approximately one and one half 
hours during which the data were taken. 

Young's theory [Ref. 11] has been verified in 
Mems. 12 and 13. This, however, applies to large apertures 
and stellar scintillation observations and would not be 
expected to apply to horizontal laser beam propagation with 
Meeciver diameters less than the first Fresnel zone size. 

Kerr has reported some letmnteed apevture averacing 
experiments [Ref. 14]. The data consisted of comparisons 
fea o2-CM receiver aperture and a 3 mm aperture at 4880 


Angstroms. Kerr points out that the predicted dependence 


ee , -X : 
Seetne 1rradiance Variance is D fmomnthe sarameter much 


ii, 





haager than oaye® and that it was not observed. Different 
puionseasSien Values of 2, 7/3 or 3 to X. The ratio of log- 
mibetude vyandances Was plotted versus strength of turbulence 
fmaoptained £rom Micrometeorological measurements. The ratio 
Seowed an increase from an average of 0.025 to about 0.080 

as the turbulence increased with the largest data spread in 
mme saturation region. {See Section II-C for a discussion 
of the saturation region.) In Ref. 15, Dunphy and Kerr 
presented typical ratios in comparison of a 32 cm receiver 
amaa 6 mm receiver. The ratios were from 0.35 to 0.70 for 
values of On tye OK co (Ey) Oe, he Oma (item at Los dard 
not always vary consistently with ce and appear to be 
Mmprosentative of the saturation regime. 

Kerr summed up the problem of aperture averaging by 
fyeeing, "An adequate theoretical description of receiver 
smoothing has not been given, even for nonsaturated condi- 
muons, [Ref, 14}. it should also be noted that all of 


mimes above data were collected from propagation paths over 


land, 


fee DEPENDENCE OF SCINTILLATION ON HEIGHT, THEORY 
Mie relationship between the refractive-index structure 


2 ; : 
ime CLON , Cy , and micrometeorological parameters is given 


by 


(24) e790. ame <—)- ~C 
i 


18 





ce wei TteMmMermature SLYUCTUre pParameter, P is the 
Meessure in millibars and 7 is the temperature in degrees 
Kelvin [Refs. 6 and 16]. As temperature and pressure depend 
on the heights they are measured at, Wyngaard, Izumi and 


Moitins [Rez. 16) developed a Semi-empirical theory to 


Meseribe the dependence of Cmriemone From this theory , 
oe ieee Pe imMiiccnOmtredmmly at tWO Separate levels, can be 
maemeten as a function of height. 

may also be written as 
(25) at = 3,.2Nb 1/5 


meere b 1S the rate of molecular dissipation and N is the 
merce Of molecular destruction of = e(i=ty->. Seal gers = Wee i 
mintiarity theory states that statistics of the mean and 
turbulent flow fields, when properly made dimensionless, are 
universal functions of only a stability parameter. By choos- 
mie height h, the kinematic surface stress ce abel. {elas 

fecal Mean potential temperature gradient ee make the 


Eyeaetstics dimensionless and the Richardson number Ri as 


mie Stability parameter, the following relations hold: 


(76 ) BP, (Ri) = —— and 
U;, 
(27) F,(Ri) = : 


hU, (30/3h) 


19 





Combining (26) and (27) with (25), 


Smee, Z 
(28) fF, (Ri) 
where . eee F., 
F,(Ri) B/S ° 
1 


Tie Sore wie nicicmec ee Ci@¢mkimonen must be further 


Helineated to eEStablish the relation between C é shige, igee 


il 
ieomadoetnaceiynpaard, et al. {Ref. 16] chose the surface 
mop hei be Lux y@ in place of ae and the dimensionless 
papame ter VimpatecmOmenimasmtie Stability parameter. 


feero detined by 
29) L=——,— 


miere Ko is Yon Karmen’s conStant and is approximately 
eoual to 0.4. 

During the day over land, when heat transfer is generally 
upward, Q is positive and . ince cteave, These are unstable 
memditions, 1.€., the buoyant air near the surface tends to 
eS\e: At night OVet inde sis cenerally negative. Over 
the ocean surface Q may be negative, zero or positive 
depending on the relative temperature influence of the 
ocean. When Q = 0, ~ = (ene se are neutral conditions which 
Mey be approached under very cloudy conditions. 


(Moet weequation (20), first convert 2% and Rain 


Taio wOtetne new variables, h, U., Q and L. 


20 





(30) hU,, 96 h 90 h 


Q oh” 1, on ~ 21‘? 
and 
(31) Ri = g,(Q) 
2°L 
Q : h 
miere 1,= Sim emt usmanrc fUNCEIOnNS OF yp only and are 
* 


Mayen in Ret. 17. When, equation (28) becomes, 


(32) Cee Tn 2/9 gH 

ime function gx Se eo aeme by, 

( Sep, g. = Asp Setifalier ie] pant for 0 > : 

and 

(34) g, = 4.9 [1+ 2.75 ()] foro <p. 
Fomeunstable Conditions, with 0 > rand - = 2 lara oo) 

becomes g, ~ xz (- _— amd (52)eas then 

5) Oa a Pia elds 


Wem stable conditions, C 
2) 5 


T decreases with height more 


Slowly than h Pormenpemeralecondit10ons, the dependence 


approaches 2/3, 


When comparing simultaneous measurements at two different 


momento in the Surface layer the following relationship holds, 


Zi 








2 -4/3 2 
Cp (1) | _ Cy ae 


(36) —7— = 
Ci 2) (CZ) Cy? 2) 


Eereeoptically determined data then, 


2 2 4/3 
Co (1) ; Cy (1) _ — 





(37) Dy a ee a 
Coy 2) Cy (2) | h (1) 
mmeereunstable conditions. If conditions are neutral 
z Z iz, 2/3 
(38) oe) ; C1) ; Cel) : Fen 
es an ee ee 
Ci (-2,) Cy (2) Con (2) Hie) 


igesaweners on this theory have pointed out certain 
Mmerprict1ons tO it {Ref£. 16]. First of all, this is a sur- 
feee layer model requiring local isotropy and may not apply 
Gieecteater heights in the atmosphere. Secondly, the relation- 
emaps used hold for long-term, or ensemble, averages of 
Meroulence properties. JTime averages are only approxima- 
tions Mom onacmolLewavetaces we Inandly, the theory is limited 
meemorizontally homogeneous conditions. This means that a 
Pevem portion of surface, with area large compared to that 
mcm typical roughness element, looks like any other portion. 
Mitoeis tO Ensure that the flow structure is in equilibrium 


mek the surface, 


ee 





C. SATURATION REGION 
io Veceri pelo 


Pe ismiieratlysdeceentca im tie Literature that 


eauwation (7), 


_ Teoh) 6 2 
Con =J a2 ae kk Z Cy 


fees not describe the retractive index structure constant 


for values of 


fae oes 2 
Coq = <CinA - nA) > > 0.8 


ne ee eee iio icles ihne effect of saturation 


dc tewmiamed Optically with e deter- 


, Zi 
1s seen in comparing Co A 


A 
mimed using values of Ge PsmeObheorncdmeron micrometeorological 


maeenods. <A one to one relation is followed until oe (op eae) 


meaches about 0.8 where it levels off and then decreases for 

increasing values of Can Peremiaigmey Reis, 14, 20 and 21}. 
Tatarski has investigated the effect of using the 

method of smooth Dereurbati ons in deriving equation (7) 

concluding that Cpe Swievaswanlimiting factor |Ref. 2]. 
Using a dimensional argument, Kerr deduces the 


jmeration 


2 ; 
Co, (Maximum) = 0.6 


Maepenaent Of wavelength, for the onset of saturation [Ref. 


14). 
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Zoe . ANE Oney, 

DeWolf has developed a theoretical description of 
the behavior of Os oe in SOndoine, he points out 
Biccethe Rytov approximation of smooth perturbations is 
meaely restricted byethe physical condition that ray bending 
Mees not derlect rays as much as one Fresnel zone off the 
Beas . 

The region of turbulence where od “ies considered 
£0 generally obey equation (7), with the best application 


Being for ann ee Nemo nell duroneG;eeglon Can be defined by 


a2 
1/3, 0% 2 ce 2 


I< Ki N 


where K is the spatial wave number corresponding to the 

mamer scale of turbulence, k is the optical wave number, 

gets the path length and Ly Peeve vouter Scale of the inertial 
Suorange. In this region the irradiance remains log-normal. 
ihe dependence of the log-amplitude variance <dx?> is calcu- 


lated to be 


(39) vere nek > 2° @.7) 1/9. 
m N 
imeecerms of an CrOmeaespilenica | wave) 
De ee t=) 1 56 2 ,-1/6 
40 = as 
(40) <a X= > 0.33 (K ie (CoA) ; 


The derivation of (39) was based on plane wave theory. In 
order to solve equation (39) So must be determined experi- 


Meteatiy, which is difficult to do. 
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Orieceauenhonusenave dttempeed to model this region, 
femeeecardace, none have suificicntly satisfied experimental 
Observations [Refs. 11, 13, and 22 - 27). 


Se GC 


on Cx 


NT 

PemIusemOrntnle shack OL am adeauate theory to describe 
mieesaturation region it 1S necessary to correlate measure- 
ments of Cu aseaeteimaimed Optically, Cue and Cy aS Gleqceiec 
fened by micrometeorological methods, Cup CinewOimmenc 
moe range goals of this project is to be able to predict, 
Om a large scale, the behavior of a laser beam on the basis 
Seemie teorolopgical measurements. The mathematical relation- 


Ships for the unsaturated region over land are given by 


ee@m@ations (10) and (24). 
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ell. ek ER IMENTAL PROCEDURE 


foe GENERAL EQUIPMENT DESCRIPTION 
1. Detector Package 

A detection system was designed and built to function 
under adverse environmental conditions and be adaptable to a 
Memtety Of experimental parameters (Figure 2). 

The detectors used were standard capacity type PIN-10 
feeacon Schottky mioroniodes Manthactumeda by Unaced Detector 
Mmeennolosy, Inc. Ihese detectors provide a linear response 


ie watts to oe Wott siaverdn NeEsrP. (KHz - .85y - 


fiom 10° 
imiz BW) of ee watts, a rise time of 10 nanoseconds and 
Meeark Current of 0.5 microaimperes. The units were operated 
mmecne photovoltaic mode with zero bias voltage and a large 
load resistance (typically 500,000 ohms). The active surface 
area (1.25 cm?) was found to be non-uniform, somewhat 
memaniscent Of photomultiplier tubes. A typical plot of 
Mmeasrtive VOltage output versus position of a focused 6328 A 
laser beam is shown in IP sen eueys. Gy, 

A 3L cm focal length lens was used in conjunction 
Met a2 thin plate aperture which limited the field-of-view, 
oes cil ticld lens to focus the input aperture on the 
photodiode. The field lens was necessary to remove a source 
Meepossible error due to either beam wander effects or 
Mmeacking errors. the lens eliminated any motion of the 


image on the detection surface and thereby ensured a 


Z0 





@emstent responsivity from the detector. The field limiting 
mete aperture was placed at the focal point of the objective 
moo Ihe aperture was removable so that its size could 
Pemchancged to meet the particular requirements of the experi- 
ment undertaken. 

bieced inmediatelyein frontesef the objiective lens 
fea 10 A filter centened at 6328 A with a 65 percent peak 
Mmeansmission coefficient. The narrow band filter, along with 
me field limiting aperture helped reduce the effect of back- 
@eound radiation. The filter mounting limited the maximum 
Mmeeulvye receiver aperture to a 23 mm diameter. In front of 
Meemriiter was an 18 leaf iris capable of being adjusted 
from 42mm to 3mm. 

ihiicwmnain teature Oofbethesdetection system was its 
mewrod Of background compensation. From Figure 2 it can be 
ween) that the detection package was made up of two separate 
M-ecilving paths. The laser signal was received through one 
path and the background radiation capable Of ibeings tiens— 
mitted through the narrow band filter was received through 
the other. The field limiting aperture in the background 
memsor was horizontally offset from the optical axis so 
mie it would not see the laser signal. A typical offset 
was such that the areas seen by the two receivers were 
aoparated by 1-1/2 meters at 4 kilometers. The placement 
Beeeenese apertures was tested in the, laboratory over a 
meen length of 280 meters prior to wse in an experiment. 


The electronics system was then given simultaneous readings 
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Pietmesdaser signal and the background radiation from a 
mierda of View as nearly equivalent to that of the laser 
receiver as possible. 
Peewee lectroni1e Components 

heplocktudiargi~an ormenie basic electronic components 
med is depicted in Figure 4. The voltage signals from the 
Meecetor package were sent to a low noise differential pre- 
mmplatier which subtracted the background radiation from 
mae laser radiation, amplified the signal and sent it through 
memmretit box to the log-amplifier. At the circuit box the 
amplified signal was monitored using a microvoltmeter, and 
meealibration signal could be fed into the system when de- 
ered, A bias circuit was’also built into the signal line. 
mec nopping the signal at the laser every 60 seconds, for a 
Sumerecond period, fthe®operator could see the amplified dit- 
Bee iicceDpctiveen the two detector signals and set the differ- 
meen ecgual fo zero, It should be noted that during daylight 
meurs it was necessary to make adjustments to the bias con- 
trol. On many occasions, the difference in background could 
memonscerved to drift, that is, the background radiation 
@manged over short periods of time. This fact was the 
mitral reason for using a simultaneous background compen- 
seron system as opposed to farst recording the laser signal 
meea period of time and then recording the background, 
memlowed by a deconvolution of the two sets of data. 

The log-amplifier was capable of accepting signals 


fee Millivolt to 10 volts. The output was proportional to 


28 





mies loscarithm (basewten) of the signal so that input signals 
oo mv, 100 mv and 10 v corresponded to outputs of 200 m, 
meoemy, and 800 mv, respectively. This output was then sent 
ema wide band power amplifier where it was amplified by a 
fefetor of ten. From the power amplifier the signal was sent 
meman Oscilloscope, where it was monitored, and then on to 
meee. M. conversion instrumentation tape recorder where it 
has stored. 

The final data processing was done back in the 
laboratory where the tape recorded signal was sent to a 400 
Meennel pulse height analyzer. Half of the analyzer storage 
Meee used to store the applicable calibration signal while the 
@ener half was used to accumulate and display the distribu- 
mon Of the scintillation signal. The sampling frequency 
feeeescet tO 16 kHz by an external trigger pulse generator. 

Uiiecmeicereuccamror the aperture averaging experiment 
moerOor the lower level of the two-height experiment was a 
me-Ne Metrologic Model TT UepmoOuidine Ameer fective Output 
mewer Of 2,3 milliwatts at 6328 A. The laser was found to 
be Mabie after a half hour warm-up period. A 6mm focal 
moeth microscope lens system mounted on the laser output 
mpenture imaged the beam on an 4cm focal length converging 
Mens, ihis converging lens was adjustable in horizontal 
position allowing the beam to be diverged a focused at 
foeee. this provided an effective transmitter aperture of 
6 mm. The laser with its associated lens systems was then 


memmtced On a four inch channel. This channel, in turn, 
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weaeemountead on a base constructed of 3/4 inch aluminum 
plate in such a manner that an operator could make fine 
fagpustments in train and elevation. 

An additional laser was necessary for the two- 
memehnt experiment. This laser was He-Ne, 6328 A, model 
S-100 manufactured by C-W Radiation and provided 0.8 milli- 
Watts power. It was fitted with an adjustable beam spreader 


and was mounted on a heavy duty camera mount. 


fee APERTURE AVERAGING 
ime Propagation Path 

iiemprepecatmonepath was 4,05 kilometers in length, 
eemepercent of which was over the ocean surface at a height of 
meoometers (Figure 5). The laser was located on the ground 
mel Of the Monterey Holiday Inn about three meters inland 
meom a 2.5 meter high vertical sea wall. The beach sloped 
mamely trom the sea wall] about 40 meters to the ocean. The 
[eeeivers were located in close proximity to a three-story 
miiding at Hopkins Marine Station. Here, the land sloped 
memcly up from the ocean about 30 meters to the receiver 
Bate. 

Data were taken on two days, 22 October and 
memcctobecr, 1973. The propagation path generally was perpen- 
Bteular to the ocean wave flow. On these two days the average 
meave neight was about one foot with three foot swells. Both 
mys were sunny with no cloud influence. The first day had 
Strong winds (20 knots) from the Southwest while the second 


day was calm with variable light winds (0 - 5 knots). 
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2 Pata We cumulation 

An earlier trial experiment across the same path 
Mmacdicated the need for multiple aperture data. The four 
Meertures used did net provide a readily defineable relation- 
ship so it was decided to augment the small end of the 
mewection with two additional apertures. Pieces of aluminum 
meock were machined to the outside diameters of 15, 9, 6.3, 
meomand 5 mm with the sixth aperture being the 23 mm effec- 
tive maximum aperture of the narrow band filters. All 
Meertures were less than the first Fresnel zone, 5.08 cm. 
fee change from one size aperture to another, the desired size 
M@eerture gauge was placed in the opening of the irises, 


Mee arises closed on them and the cauge removed. The move- 


& 
emt Of the arases was tight enough to ensure that there 
meaoenro Chace im the opening once the size was set. 

Prem or prem ren discovered in the trial experiment 
fas that S€parate aperture data accumulation periods could 
mon De Correlated without the possibility of significant 


Perors due to the change in C TOeciIEecunvent this prow rem 


N° 
Beother detection package was required to provide a continu- 
ous measurement of Cy, eaminxed anerture. lhe detector 
package described in Ref. 6 was modified to match the 
ieximum aperture size of the variable aperture unit and a 
Suatable field limiting aperture plate was machined to match 


mae Other unit. The two packages were tested in the corri- 


men laboratory to see if they would yield identical oe 
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Mamie etne 25 mm dpertures. The systems agreed to £1.53 
percent under weak turbulence conditions by placing them 
peecre by Side. 

Micmac cnmone22 anc 27 October, then, consisted 
mera continuous monitoring of Cy With a 23 mm aperture while 
Mee ine tne aperture Size on the other detector package, Each 
eomplete data sequence started and ended with a run at 23 mm 
on the variable aperture package, taking about 14 minutes. 
Maen yun Was of two minute duration with the laser signal 
being chopped every minute for background bias adjustment. 
Mieesionals trom the receivers ran through parallel electronic 
menponemnts to the F.M: tape recorder where they were simulta- 
meeusiy taped on different channels. An additional voice 
em@amnel was utilized to make pertinent observations, in- 
Sading time calibrations for each run. Six data sequences 
mere taken on 22 October and seven on the 27th of October. 

De IDE se\ erahb here ckonol 

Pie wmiimie mh HOmltOl ves tae ntaped Signals were fed 
mic<o an OSCilloscope for monitoring purposes and from there 
mico the pulse Hemutematiaiy cere PEP 400. ihe first signals 
mecorded during each experiment were the calibration signals. 
miese Signals were 50, 100 and 200 millivolts at the input 
Side of the log-amplifiers. They were stored on 200 channels 
at the PIP 400 in such a way that they could be superimposed 
over HicmeoOc mieenoity Gdistributions aiter each run was 


mecumulated on the other 200 channels. 
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iewimrecicttysGatstributions were log normal, as 


Mecaictea by theory [Ret, 26]. The desired quantity was 


Z 
Nol 


@e@eantity displayed on the PIP 400 was 


C ~ TING WelIeLeeS Oi epost ‘Wes esi pu tlon. Une 


1) alog,,! = 0.434 alnlI = blnil 


mrapeond and b are GonsStants. The full width at half maximum 
of a normal distribution is Xea/2 = oe NeTe Bene , 


(42) g = 0.425 Xe 79 


Both o and Xe 72 are measured in units of whatever 


Mmemaistributed, in this application blnI. Now 


2 


I Z eeZ laps 
(43) <(bln zs) eee OMS) 2->~b'C 


at 





memeectromim tie right units, a scaling factor, 5S, based on 


mae known calibration signals was used. 


b{inI,-1nT, | 


a ie o> erm 
. X, xy cm 


Then the variance of the Greermputrom could be found by 


measuring Xe1/2 in cm and using 


Bees _ U2 bin(15/1))X eq 2 Cem) 


X, - Ay cm 


mach in turn yields 


5 





b [0.425 In(1,/1,)Xpy 7] | 
——e es Bei 


ond finally, 


; 2 
OA 5 In(14/T))X¢q 2 Cem) 


QT X, - Xy cm 


(46) c 


EAcheacdiSstniDielom, Withee the Calibration signals 
Superimposed, was photographed. Xe1/2 was then measured off 
mame photograph using a vernier calipers. Maximum error 
mmeroduced by this technique was +2.7 percent. Cy was 
computed using equations (10) and (46). 


fe LVO LEVEL EXPERIMENT 
Meeeconoimeaieteorological®and Optical Effort 

MicmpumeOse Olsthus Cxperiment was to test the theory 
Meeecitcd 11) Section 11.8. over the ocean surface. The same 
propagation path as described above was utilized for the 
lower level. The'2.3 mw laser and its corresponding re- 
ceiver were 4.9 meters above the ocean surface. A second 
System using the 0.8 mw laser measured scintillation at a 
Mereht of 12.3 meters, directly above the lower path. 

The Research Vessel ACANIA was fitted with the 
ferecorological instrumentation and anchored about 15 meters 
south of the mid-point of the propagation paths. Data were 
meen trom 1945 to 2221 (PDST) on 18 October. The wind 


was from the South at 3 - 5 knots; waves were about two 
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feet hich withetwo foot swells. For most of the experiment 
the area was covered with low clouds, 
2. Data Accumulation and Reduction 
ae Meteorol oc eal 

UiicnomPleiwotmunc Cxperiment was performed by 
members of the Meteorological Group. Meteorological instru- 
ments and measurement techniques were essentially the same 
Pemeerose reported im Ref, 29. An instrumentation mast was 
gnstalled on the bow of R/V ACANIA. Two amplifier bridge 
SameuitsS, (Sylvania GIE Thermosondes), each with two 60 
om platinum ware sensors were mounted at 4.525 meters and 
fee. meters above the water line. C., eould then be obtained 
Pyeemicasuring changes in the resistance difference between 
mre two Wire sensors spaced 10 cm apart. The resistance 
mmerercncesis proportironal to the temperature difference. 
Peyoltage signal proportional to AT (rms) was recorded on a 
Strip chart as a function of time. 

In processing the data, raw values of AT (rms) 
weme abstracted from the recorded chart. These values were 
mem corrected for calibration errors in the nominal resis- 
mamec Of the platinum wires, 


Cr was then calculated using 


Cor = true AT ai eee 


mere ris the probe separation of 10 cm. From this, Cy 


was found using equation (24). 
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De Mi meiitinoms were CStimated from observa- 
imioicomOtmeremeradicnt formed by the sea water temperature 
and air temperature sampled continuously throughout the 
Peeperiment. Quartz thermometers were placed at the lower 
level and on the ocean surface. 

bee Opt 1.e@agl 

(Cee Menpiemeetised CO, receive and process the 
data was the same as for aperture averaging with one addi- 
Oren. Alt the electronic components were located at the 
Meier Teceiyer Site to limat the number of operators to 
eee liis required the use of a long coaxial cable from 
mae Upper receiver to its differential amplifier. To ensure 
mem losses were encountered, a cable driving cathode-follower 
Merectit was built usine field effect transistors. The 
Sim cult performed the function of subtracting the background 
Signal from the laser Seitiidomnecl beds driving the long 
Gable. 

Tiemireidmlinicimngeaperture plates were machined 
Mmmeover dad J.l5 meter diameter circle at the 4.05 Km distance, 
This, together with 3.7 meter beam diameter sizes at the 
meceilvyer ensured that each receiver was only seeing the out- 
Bem OL 1tS respective laser. 

VitoonvetemaGgecumilarted £0r an average of three 
maces per run. Ali data reduction was accomplished in 
the laboratory using the same techniques as described in 


mection [II.B.3. 
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ie EPREQUENCY SPECTRUM ANALYSIS 

Representative samples of data were chosen to compare 
mie erelataive power spectra of scintillation. These samples 
include upper and lower level data from the two-height 
experiment, comparison of the two detector systems when side 
Myeside (12 cm apart), and a complete aperture averaging 
Sequence. 

The recorded scintillation signal was sent through a 
wave analyzer with a band-width of 10 Hz. The lower 
meequency Jimit of the analyzer was 20 Hz. Each signal was 
m@imecred at 10 Hz intervals from 20 to 100 Hz, at 100 Hz 
meamervals irom 100 to 1000 Hz, and at 1500, 2000, 3000, 
4000 and 5000 Iiz. The output was sent to an X-Y recorder. 
mie power was determined by integrating the output curves 


with a planimeter. 
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IV. EXPERIMENTAL RESULTS 


fee «=APERTURE AVERAGING 
1. Preliminary Investigation 

[fic iiieanvestmemeton Of the E€flect Of aperture 
averaging on or using diameters smaller than the first 
Fresnel zone was conducted in the laboratory. The C-W 
Radiation 0.8 mw 6328 A laser beam, slightly diverged, was 
meme down a long corridor and reflected back to the receiver 
ivea plane mirror. Total path length was 286 m. Turbulence 
meas provided by seven large forced draft ventilators which, 
admittedly, was not representative of turbulent flow in 
me open AGIOS NOSES Diemilevetemetess. did produce a log- 
mmeensity distribution that proved to be normal. 

Each sequence of measurements contained data taken 
mene apertures Of 3.5, 9, 18 and 23 mm. Relative values of 


C Meme catculated by using Cy obtained with the 3.5 mm 


LA A 
werrire 2s the normalizing factor. The results of five 
@oea Sequences and a composite plot are shown in Figure 6. 
Each data point is an average value taken from five or more 
mesecond sampling periods. The degree of aperture averaging 
mound here indicated that further investigation on a full 
peale path was necessary. 

Peete rey Day Experiments 
rer eOeroDer 1975 
Meowereovlous.y mentioned, this day was somewhat 

unusual in that the winds were quite strong and were from the 


penin, 1.6., over land. 
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ThewGem@EnOhenecervyer SyStem, operated with a 
fixed 23 mm aperture, was used to continuously record data 
feome 1555 to 1751 (PDST). Figure 7 is a plot of the values 
of Con Sotamied Whenmuence Vandablewaperture system was set 
mem25 Mm VS. time. Only one data run provided identical 
values of Con EVOMMEGUNneceeaVver Systems. For the remainder 
@eethe runs, Con Witecmormtdinedewiti the variable aperture 
Eeeotem were as much as 12 percent lower than the control 
Byotem Values. The average deviation was on the order of 
mepercent lower. Some factors that may be responsible for 
mis ditherence are beam wander, instrumental error and 
meemtion Of the recenavers with respect to the beam axis. 

Assuming that the calibration set in the 
laboratory was still valid leads to an investigation of 
@ecombination of beam wander and position relative to the 
beam axis. (Beam wander was measured at approximately 15 
Manmute intervals by Beall and is reported in Ref. 9.) 
re intensity distribution of the beam in the absence of 
may turbulence can be conmsidered to be Gaussian. If a 
Petector is placed in the beam and the beam is allowed to 
Meander randomly, the intensity differences sensed will 
mepend on the magnitude of the beam wander. If the detector 
Mmempon tie Original axis, the highest intensity sensed will 
mamigen the beam Genter 1S coincident with the detector and 
the lowest intensity will be when the beam is at its maximum 


meemrsion trom the Omweinal axis. “lf, however, the detector 


59 





iemplaced atea point turther off the original axis than the 
maximum beam wander, the maximum intensity sensed will be 
at one limit of the maximum beam wander and the minimum will 
merat the Other. The intensity difference, then, will be 
mien larger than that sensed when the detector is on the 
beam axis. ‘ion Duopagatime 1m turbulent atmosphere, 
these intensity differences are coupled with those caused 
by the medium inhomogeneities and increase the apparent 
weuntillation, 

Onm@mitmemecday the position of the receivers in 
ere beam was found to be about 0.35 m off center. The beam 
diameter was 2.74 m and the distance between sensors was 
erie ihe Control (1.e., reference) system was further 
away from the beam axis and, as shown in Figure 7, always 
meee tne Jarger value of CoA: 

H@m@ctes sO nas Shown theoretically that the 
fmerlance Of the log-amplitude fluctuations, Cons increases 
with radial distance from the axis of a laser beam when the 
Mmeerture 1S located in the near field. The 4.05 km propaga- 
Bien path Mime mer adtismteramomnatler aperture placed the 
Meemelvers im the far field, A similar effect to that 
described by Ho was, however, observed in the corridor labor- 
meory under far field conditions. This StnoGr me anplrtved 
by beam wander, appears to explain adel eeerolly fen wre: dick — 
ference in Con values Obtained . 

MicmpurpOsc Olmene Control System was to 


continuously monitor the changes in the atmosphere so that 
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the values of CoA obtained with the smaller apertures could 
be compared to a known reference. Because identical values 
of C,, were not obtained Micietm@e 25 Mn apertures, at the 
beginning and end of each sequence, a linear interpolation 
fas Carried out to obtain BDaseline data at intermediate 
memes, Ihis interpolation process consisted of comparing 


the difference in C vanes observed at the beginning of 


RA 
me Ucnece COntic@cittorenee 11) Co, Values OM stenevicw ra ee aaie 
pecumo. tne sequence, Define these values as AC), (ty) and 


AC, ,(t.), mdwase tire Gare of change of AC, ,(t). Then 


: AC) y(t ) - AC, , (ty) 


ie - ash 


ire interpolated value of Co, (25mm) for the variable aperture 


meeweem at a time t in the interval t, < t < t. would be 


il Z 


eeual to 
Co, ft) + [AC , (ty) = F(t-t,)] 


where Co, (t) is the value of CoA obtained with the control 
fs tcm at time t. 

Piiieeswomanad J depict tlhe relationships between 
Co, WD) as seen by the respective apertures and Co, (25 mm) as 
Seecrmined by the interpolation process. These figures 
Show that Co, (D) does not consistently increase and decrease 
with Cy, (23mm) apa st tatetnc value of Cy, D) - Cy, (25mm) is 


imeeier CONnSiStent for a given aperture. 
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Pit omammlOoteor the diiterences between 


LA 


meriterpolated from the 23 mm aperture measurements. The 


mne values of Co 6D) for the respective apertures and C 


murye rises with decreasing aperture size until 4.8 mm 1s 
meeched. The average value for the 3 mm aperture was slightly 
Hess than for the 4.8 mm aperture. 
bee 27 BOGEODe. 2197 5 

This day was sunny with very light and variable 
meoas, the same procedures were followed as above. Figure 11 
mea plot of the 23 mm values of Con as determined by both 
meses, — Once again, the two systems did not yield identical 
memes. Lhe control receiver system was closer to the average 
Mesettion of the beam axis. The difference in the plots may be 
Gealitatively explained using the same arguments as in the 
Meevious Section and adding the fact that the overall value 


of C,, was larger. (Beam wander measurements were not taken 


N 
emethis day. ) 

bome Pnterestine visual observations of beam 
memaer were made. The laser was first sighted in on the 
receiver oreerneounaehz2u0 (PNST) but kept driiting vertically 
downward. Because of the laser elevation adjustment mechanism, 
myeevertical drift due to mechanical slippage would have 
Caused the beam to move upward; therefore, the beam motion 
Peeedue to the atmosphere. The rate of drop was about one 
fewer per 10 minutes. The beam stabilized from about 1330 


Semeeoo0 {PPSl) at which timewit started to drift upward at 


meeut 0,s meters per 10 minutes. A possible explanation 
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Moldepe large seale density changes, or more likely density 
gradients. The sunny calm conditions were conducive to the 
formation of a2 non-disSipating warm air mass over the cooler 
ocean surface. This could form a density gradient with the 
fiemcity greater nearer the suxface. The beam would then 
mena tO bend downward. As the air above the surface cooled 
down later in the day, the beam would move in the opposite 
aerection. 

Figures 12 and 13 show the relationship between 
the aperture averaged system and the extrapolated 23 mm 
moves. Figure 14 15 a plot of the difference between the 
extrapolated 23 mm values and the values obtained using the 
mawerous apertures versus aperture size. The shape of the 
meve 15 Very similar to the one obtained from the data of 
miee22nd {Figure 10). The peak of the curve in Figure 14 
Mmoeecz.0 times as large as the curve in Figure 10 indicating 
that the aperture averaging effect for apertures smaller 
mien the first Fresnel zone may be functionally dependent 
on Cy. 

- ¢. Composite Data 

Paces LP onieeciov that the lareest relative 
Malues of Co, Were obtained when the 4.8 mm aperture was 
used, Based on this, it was assumed that the 4.8 mm aperture 
esaproximated a point receiver better than the other apertures 
used. 

Hieouccre torrid ag relationship between aperture 


averaging and Co: HmiimemNOleattOme process, Similar to that 
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carried out on the 23 mm apertures, was undertaken to 
establish a continuous set of values for C,, as seen Dye 
imo mm aperture, the ratzos of Co, CD) /C,, (4.8 mm) were 
mmtcilated, averagzed fOr each collection aperture each day, 
MoeplhouLcd a5 aueunctHen Of aperture diameter in Figure 15. 
Values of Co, (D/C, , (1mm) obtained from curves published by 
memed, et al. {[Ref. 7] and reproduced as Figure 1 are 

Meer ted +01 Comparison. Linear regression lines were calcu- 
[iereanand plotted £0r all three data sets. The slope of the 
Meeression bine for the 27 October data was steeper than 
Biecewcorresponding to 22 October and Fried's data. The 
imenes Cross Im the vicinity of the point Con CDI/Cy y equal 
monc.95 and aperture diametcr equal to 10.6 mm. Figure 16 


moe plot of the absolute values of the slope + vs Cy, as 


A 
mem by the 23 mm apertures for the regression lines plotted 
mmeereire 15, The line was calculated using the least square 
method. (Fried's data were adjusted to apply to the 4.05 km 
eee) Using this limited amount of data a point collection 
aperture jgeMlley ene Con Sertideoc westimatcdsimethe L£ollowing 
Hemmer, Entering Figure 16 with the observed value of 

Co, 023 Mie telassticeinyerse Shope to be uSed in Figure 15. 
meplying this slope in such a way as to cross through the 
eaperveally observed point of intersection of the regression 


mimes will give a value for the ratio Co, (23 mm) /C selon 


Nes 


miecturn yrelds Cone 
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few iosODwlOuSmenat additLondl data are needed 
for LOormudate a reliable mathematical relationship between 
Peerture diameter and CoA or Cy 
da Peeeteney, spectrum 
Chi teovdmrcewoiinas calculated the theoretical 
temporal frequency spectrum, W(f), for a spherical wave. He 


solved the plane and spherical wave asymptotic forms in 


ey: 
Berms of CoA? i and 2 


where £, = v/(2raL)*/?, 


at Mec ti ime HOG perpendicular te the propaga- 
tlon Daath. and 


Q = tf ee 


For a plane wave 


ae ee 4/3 
(47) NCE) = 0.44 (Cyd GI Ct0.2707°), 2 << J 
(48) W(£) = 1.14(Cy 4) 675) etd eae 
O 


Mmieene Spherical wave case 


a 2 1 4/3 
(49) WG Deed. S (Ce Fee? aan? pees 
(50) W(f) = 2.81 ae ii  R s3 2 
O 


eiese solutions are valid for a vanishingly small aperture. 
iieagskiv (rer. 2) calculated the effect of 


meecivyer aperture averaging on the frequency spectrum for a 
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plane wave model. This effect shows up in an additional 
feoemultiplied on the kérnal of the integral, which corre- 
Ponds tO agpartial Suppression of the normalized rluctua- 
muons. this factor is 


5 Z 
Tree meant /V ) 


3. 


fel) a 


Rae eet / V 


where f is the frequency, v is the wind velocity perpendicular 
moethe propagation path and K is the spatial wave number. For 
mee oevy ZR, tnis term 1s small for all K, For high frequen- 


mes the aperture averaging weakens the effect of inhomogene- 


metecson all scales. If £°<< v/27R, the factor is close to 
Party for k << 3 and becomes small for K >> 5 » AUN Clennes= 


bution of large-scale inhomogeneities in the low-frequency 
meme ol the Spectrum 1S virtually unaffected, whereas all 
Scales which are small compared to the aperture radius are 
suppressed, 

i hemicepoOWcr Values bor cae aperture were 
@orained by dividing the power found at a particular fre- 
quency by that found at 20 Hz. From equations (49) and (50) 
meecan be seen that the ratio W(L,)/W(L,) for a vanishingly 
eiebt aperture is a function of £ {> and f, only )eihre 
Meerture dependence is contained as a multiplicative term 
mad tirst order Bessel function as shown in (51). The 


Mependence should not be lost in taking the ratio 


WCE, J /W(£,); therefore, when comparing relative power vs. 
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Peienea setems, a diiterence Should be detectable due to an 
Meexture change. If the data samples were not taken simuita- 
Mmeously, Changes in v could also contribute to any differences 
epserved, EStimates were made of v because the propagation 
met could not be properly instrumented. 

Ome 22 50Gbober, vy Was eStimated to be 16 knots 


(8.24 meters/second). Defining 


(52) £. = 2m 


mren for this day - =e ll) on (gil) ae ou re £ (6.3 mm) = 417 Hz 
and i Coe —- iideazee Freure 17 1S a comparison plot of 
relative power vs. frequency for data obtained with a 23 mm 
aperture and for data simultaneously obtained with a 6.3 mm 
eperture. 

it VUMiZ wire brelativie power Curves in 
Peete jf are converging. For most frequencies sampled, the 
fem aperture receiver showed a smaller contribution to 
melative power. This is in agreement with the aperture 
dependence and with the previous statement that all scales 
meinhomogenecities which are small compared to the aperture 
radius are suppressed. 

On 27 October the winds were light and variable. 
Tf an average value of 1 knot is assumed for v, a Sree 2 
fee Rk) and £, = 4 Hz. The value of ie einem re Oils oop 2 
ome Oieaperture to 7 Wz for a 235 mm aperture. Figures 


meeand 19 are plots of the relative power vs. frequency 
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PommieneeScqucmicGe OF apertures. Ihe plots are similar in 
eiape tO Figure 17. 

Pitinemhomiae plotted Lrom data received by 23 mn, 
pemm and 4.8 mm apertures. The values of i hoOLev edual tO 
Mero are / tiz, 18 Hz and 34 Hz respectively. For f£ < 90 Hz, 
Variation with aperture size of the relative power in any 
mseequency band 1S unpredictable and is probably due to slight 
feeracions in v for the different sample periods. Sample 
periods for the different apertures were four minutes apart. 
memes - 90 Hz, the smaller apertures show a larger contribu- 
prem tO relative power and the values tend to converge on each 
miner @5 they approach zero at about 500 Hz. This is in 
mereement with Tatarski's model for f >> f The temporal 
Weariations in v would have less effect in this region. 

aie Lomdianed Wene received by 15 mm, 06.3 mm 
m@ees Mmm apertures. The corresponding values of ft it ae 
Peer cOsleknotware WhHz2,.26 Hz and 55 Hz. The same be- 
Mmavior Of the relative power vs. frequency is observed as 
maeeicure 1§ cxcept that the region where ordering of 
W(£)/W(20Hz) is seen appears to be f > 300 Hz. It would be 
expected that this region would begin at a higher frequency 
mento the requirement that f >> f ANCdetivee values. Of f. 
Zoetnis Praphi are generally larger than for the apertures 
Meead in Figure 18. | 

DiwuUnee7 io? CONMparison Of relative power vs. 
imeedqucimey as received by the 6.3 mm aperture on 22 October 


ameeez/ October. Ihe values corresponding to 22 October are 
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Pecimecmuc consistently larger than those of 27 October. 
iivseaititerence is due to the significantly larger value of 
Parone 22) Octoper,. A comparison of this figure with Figures 
memonGe lI shows that the effect of v on relative power is 
Peeater than that of aperture size within the limits of 
these experiments. 

Figure 21 is a plot of relative power vs. 
Bmequency as seen by the 23 mm apertures of the variable 
Mt weeiLewoyotcm and ene fixed aperture control system. The 
eenerol system shows slightly lower values of relative power 


mem most of the frequencies sampled. 


ieee DEPENDENCE OF SCINTILLATION ON HEIGHT 
leer rca | NMeasmrements 

Gpercal measurements were taken from 1944 to 2219 
cemeeoeOctobper 1973 (PDST). The plot of Con WE 4 Fe iviley aoe (ei 
Pwo léyvels 1S Shown in Figure 22. Both receivers utilized 
meoemmeaperture. the values obtainéd at the upper level for 
miesiirst 40 minutes are not considered to be representative 
because eit resend neainine the beam axis on the 
jmeeciver, Ihe lLamge gap in upper level data from 2023 to 
meds was due to making beam adjustments. The data from 2105 
momz206 Shows a definite relationship between the two levels. 

lables ledisomlays the average values of Con Oe eee 
mmrrerent averaging periods and their relationship to the 
Mmements of the two systems. Period (a) was from 2105 to 


mee, 1d) trom 2140 to 2206; and (c) from 2209 through 
PZ) (PDST). 
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TABLE 1 


Period. (2) Period (b) Period (c) 


C,, at h(1) = 16 ft. 0. 360 0.349 0.271 
Con cemee = 40 ft. Oe oz G2 197 e204 
CCD AMOL) Il 3.44 3.44 
(h(2)/h(1)) 2/9 1.85 
e*, (1)/C2, (2) asa 3.30 ey 
2A LA ; ° nae 


MomGEseusscad sili cetlon 11,82, during stable 


meteorological conditions Cn (and therefore C.*) deéereases 


mich height more slowly than pool Under mecutral conditions 
mime dependence approaches mn ore Pidmidchmilo tan le Conditions 
-4/3 


~ waries as ih L[Omes tobremwGOndilL1OnsS with data taken 


Pam caneously at h(1) and h(2) 





ae are itt) | ay 


(38) G7 (2) imc,“ (2) (2) on 


For Unstable conditions 





co Gb Ca n¢2))4/3 
:, oa) — oy : 
bat) Cy (2) IPyligsh,) 


Ore Cremresuits displayed ain Table 1, period (a) 


4/3 


Emows that (h(2)/h(1)) agrees with eet /C* (2) within 


Mepercenc; Lor period (b) the agreement is 4 percent. For 
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ZS ’ 2 2 
period (c) (h(2)/h(1)) / agrees with Coy C1)/Co, C2) to 
within 4 percent, as would be expected for neutral 
conditions. 
Table 2 shows the same data with aperture averaging 


@errections added using Figures 15 and 16. 


TABLE 2 


Period (a) Period (b) Pemiod™ (c) 


CoA ae 1 (1) Og Ea O2425 Gao 

Cy, at h(2) 0.226 0.226 0.240 
mm) /h(1))!- 3,44 3.44 

(h(2) /n(1)) 27° | 1.85 

1) /C~. (2) 48 % Bs 1.82 
an : . 


memuod {4) shows agreement to 9.6 percent with the predictions 
Mereunstable conditions. Period (b) is within 2.5 percent of 
agreement with unstable Eomdtt1onseand period ({c) 1S within 
i. 5 percent Preemeni wihtrmeutral Stability conditions. 
oa Gptical bie @qucher a opeetird 

Peete etSea Dlot Or relative power vs. frequency 
Peco oy tule two levels at 2154 (PDST) 18 October. The 
meow ts very much like that of Figure 21 where the two 
Meweetton systems were placed side by side, From equations 
er, (50) and (51), the theoretical temporal frequency 


eeeetrum, W({fL) 1s a function of e fe Le ene Msc sions 


one OF 
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mato W(t)/W(20HzZ} varies only with es giieae Lor a) civen 
frequency, the relative power values are dependent on the 
mercctive Maenitude Of vat the sensor level. Figure 23 
mmdicates that the value of v must have been comparable at 
the two levels on 18 October and less than v on 22 October. 

Wea Mereconuo LOvueameMedsurements 

Romg@etined iMsoeetion fF.B., when the surface 
Meamperature tlux, Q, iS positive, the air mass is unstable, 
Meeaual to zero corresponds to neutral conditions and Q 
Mmewative corresponds to stable conditions. The temperature 
difference between the sea water and the air at the lower 
iwiememeeo2> Meters above the Surface, is plotted in Figure 24, 
Conditions were unstable from 2105 to 2200 when the ship left 
M@ieeared. Although quantitative measurements of Q after 
2200 were not available, low lying coastal clouds were 
Met@eine in throughout the evening with a ceiling of 200-300 
meet by 2220. The increasing cloudiness would tend to 
mmomee Conditions toward neutrality and the corresponding 
ayo relationship. 

A malfunction in the lower level temperature probe 
System prevented the comparison of the two level optical data 
with meteorological data. The difference of 3.9 meters in 
the upper levels of the optical and meteorological sensors 


Mmeaevenuca 2 point by point comparison of C,, values. 


N 


SZ 
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V. SUMMARY AND CONCLUSIONS 


A. APERTURE AVERAGING 

Se ourlcCant apetreumc averaging was observed in an 
mIerocca laboratory fOr apertures less than the first Fresnel 
mew scizZe. lwo ELleldvexpermments, conducted under different 
turbulence conditions, indicated aperture averaging, for small 


Mier tures, 15 dependent on C Damar cOlleCetcdmpy Fried, 


Ni 
Mevers and Keister [Ref. 7] lends support to this conclusion. 
Both field experiments show that a 4.8 mm aperture can 
approximate a vanishingly small receiver and should be used 
meme vor pOosstiple. = None Of the existing theories adequately 
describe the observed data. 

iitectroneecied primitive method of predicting Co, CO) 


from Co eWserved by @ 25 mm collection aperture appeared 


A 
Meme succecesstul for the two level experiment data. Addi- 
Mm~onmealeaperture averaging data 15 required to make accurate 
poeedictions of Co (0). Uieecapanmility to make such predic- 
ions would greatly reduce the Precis tewner Unit »area 
mained to collect good scintillation data. 

The observed frequency spectrum is in general agreement 
Peace predicted trends in Clifford's model for a spherical 
meee (Ret. $1] and Tatarski's aperture averaging modifica- 
mon term aS applied to a plane wave [Ref. 2]. The effect 


feeonlallecr apertures seeing more relative power in the 


higher frequencies than larger apertures was observed. For 
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miemEiec OL apertures uscd (3mm to 23 mm) the effect of 
mpeniire Size Om relative power at a given frequency was 

fe eMeeonbe Much stess tmam the vy dependence. An increase in 
Peshows Up aS an increase in relative power in the higher 


frequencies. 


eee DEPENDENCE OF SCINFIELATION ON HEIGHT 
2 
N 


(Oneunscetanle meteorological conditions 


Wigaand kzumieanmd Collins [Ref. 16] predicted C 
4/3 


was dependent on h_ 


“ays LOmencUmnaImNckeCOnOlogrTcal Conditions in the 


and on h 
mire xt OF Long time averages, near the surtace, over land. 
ites data presented in Tables 1 and 2, though sparsc, are 

mamecxcellent agreement with these predictions and indicate 


enat the semi-empirical Pic Oils applicable to propagation 


maeens Over the ocean surface. 


fee EXPERIMENTAL PROCEDURE 

PCeUnulaeronwOrescantrilatvon data for the aperture 
meee Ine experimenes and the two level experiment pointed 
Pee Some potential sources of error. Beam wander makes a 
significant Teel oMmmOMwEO utes SCINtitlatilon observed by a 
@etector. If meteorological measurements are to be used to 
Peecict the amount of Scintillation a laser beam experiences, 
an additional term should be found to compensate for the 
effect of beam wander. 

iiemvalue of Cy ies eeOneilnyedenmendent On the height 


of the propagation path. This fact can be used to obtain 


mearec ranges of Cy On any given day. When making comparisons 
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of C and C. je eOuldebe takem to emsure all .optical 


NO VT 
Ser pOnenes ana MetcOrological sensors are on the same level 
because the -4/3 and -2/3 scaling laws are not applicable 


me short time average Comparisons, 
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